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Executive  Summary 

A  series  of  studies  were  conducted  using  a  rectangular  rocket  chamber  burning  LOX/mcthanc 
propellants  to  investigate  the  effect  of  the  temperature  of  the  methane  on  the  initiation  of 
combustion  instability.  The  rectangular  rocket  used  in  these  studies  employed  LOX  centered 
coaxial  swirl  injectors.  The  rectangular  injector  used  in  these  studied  provided  16  possible 
injector  element  locations  in  a  linear  array.  Variation  of  injector  element  location  showed  that 
injector  elements  placed  near  the  pressure  anti-node  of  the  first  transverse  acoustic  mode  of  the 
chamber  resulted  in  higher  pressure  oscillations  than  injector  elements  located  near  the  pressure 
node  of  this  mode.  The  methane  propellant  temperature  varied  over  a  temperature  range  from 
-70  to  - 1 00°F  (-56  to  -TS^C)  for  the  cooled  studies  and  from  -6  to  +40°F  (-2 1  to  +4°C)  for  the 
ambient  temperature  studies.  The  results  obtained  from  monitoring  the  pressure  oscillations 
using  high-speed  pressure  transducers  indicated  no  difference  between  the  cooled  methane  test 
cases  as  compared  to  the  ambient  temperature  test  cases. 
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1.0  Introduction 


Combustion  instability  is  well  recognized  as  a  major  risk  in  rocket  engine  development 
programs,  particularly  for  liquid  propellant  rocket  engines.  Oeeurrenees  of  combustion 
instability  events  have  been  observed  in  numerous  engine  development  programs  dating  from  the 
early  1950’s.  Particularly  troubling  is  that  no  fundamentally  based  design  methodology  for 
predicting  the  onset  and  sustenance  of  combustion  instability  in  rocket  engines  exists. 
Furthermore,  because  chamber  size  is  important  in  defining  the  acoustic  modes  of  the  chamber, 
incidences  of  unstable  behavior  are  often  not  observed  until  the  full  scale  engine  is  tested, 
usually  near  the  final  stage  of  the  program.  From  a  cost  and  schedule  perspective,  this  is  the 
worst  engine  scale  at  which  to  address  a  critical  problem  that  places  the  entire  development 
program  at  risk. 

A  major  reason  for  the  lack  of  a  fundamental  understanding  of  combustion  instability  is 
due  to  the  complexity  and  multi-variable  nature  of  the  phenomena  itself  Initiation  and 
sustenance  of  combustion  instability  has  been  shown  to  be  sensitive  to  a  number  of  potential 
driving  mechanisms  including  propellant  injection,  inter-element  injector  flow  interactions, 
atomization  to  form  drops,  secondary  atomization  of  drops,  drop  heating  and  vaporization, 
mixing  processes  involving  the  drops  and  gases,  mixture  ratio  distribution  in  space  and  time, 
chemical  reaction  and  gas  dynamics.  Different  stability  behavior  is  also  seen  for  different 
propellant  combinations  (e.g.  liquid  oxygen/hydrogen  (LOX/H2),  LOX/kerosene,  storables,  ete.) 
as  well  as  for  different  injector  geometries  (shear  coaxial,  swirl  coaxial,  impinging  jet,  pintle, 
ete.).  Thus,  it  is  not  clear  that  a  single  mechanism  is  responsible  for  the  oeeurrenee  of 
combustion  instability  in  every  rocket  engine  in  which  combustion  instability  is  reported  or 
whether  multiple  mechanisms  are  eoneurrently  present.  One  observation  is  clear,  it  is  rare  for  a 
rocket  engine  development  program  to  not  experience  rocket  instability  problems,  and  the 
current  understanding  to  address  such  problems  is  largely  empirical  in  nature.  It  is  also  clear  that 
analysis  capabilities  to  understand  the  initiation  and  sustenance  of  combustion  instability  must  be 
improved  if  high  performance  reliable  rocket  engines  are  to  be  developed  in  a  cost  effective 
manner.  Furthermore,  it  ean  be  stated  that  gaining  this  needed  capability  will  not  result  from  a 
single  research  effort  but  will  rather  only  be  achieved  by  a  systematic  series  of  studies  that 
investigate  the  underlying  fundamental  processes  that  drive  combustion  instability. 
The  overriding  objective  is  to  have  a  continuing  improved  understanding  of  combustion 
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instability  such  that  the  best  possible  tools  arc  available  when  rocket  instability  problems  arise. 
Ultimately,  one  would  like  to  have  a  eomprehensive  understanding  of  combustion  instability,  but 
history  has  shown  that  such  an  understanding  is  elusive.  Yet,  with  each  study,  new  key 
understanding  is  gained  that  contributes  to  aehieving  the  desired  goal  of  a  eomprehensive  model 
that  provides  a  design  methodology  to  avoid  oeeurrcnces  of  eombustion  instability. 

Given  the  ubiquitous  nature  of  the  oeeurrenee  of  eombustion  instability  in  liquid  propellant 
roeket  engines,  nearly  every  review  of  eritieal  areas  for  whieh  improved  understanding  is 
required,  ranks  it  at  or  near  the  top  for  needed  future  work  [1].  Historieally,  eombustion 
instability  has  been  most  problematie  for  LOX/hydroearbon  (partieularly  kerosene)  and  storable 
propellants  sueh  as  nitrogen  tetroxide  in  eombination  with  the  hydrazine  family  of  fuels. 
Engines  using  LOX/H2  propellants  have  experieneed  considerably  less  incidences  of  combustion 
instability,  but  are  not  totally  immune  to  such  problems  [2,  3].  However  the  use  of  liquid 
hydrogen  introduees  serious  systems  related  issues  due  to  its  low  density  leading  to  large  fuel 
tanks  and  its  low  vaporization  temperature  limiting  its  storability.  Reeently,  renewed  interest  in 
lower  operational  eosts,  higher  propellant  density  and  less  ehallenging  thermal  storability  has 
renewed  interest  in  methane  as  a  fuel  to  be  used  with  LOX.  Methane  offers  the  advantage  of 
being  a  “soft  cryogen”  as  its  vaporization  temperature  is  signifieantly  higher  than  that  of  liquid 
hydrogen  (-259‘’F  (1 1 1.7  K)  vs.  -423^  (20  K)).  Thus,  both  the  Air  Force  and  NASA  are  looking 
towards  LOX/methane  for  upper  stage  and  in-spaee  propulsion  applications. 

Given  this  reeent  interest  in  the  LOX/methane  propellant  eombination,  the  present  study 
foeuses  on  combustion  instability  in  LOX/methane  systems.  In  partieular,  emphasis  was  given 
to  investigating  the  effeets  of  fuel  temperature  on  the  onset  of  eombustion  instability  for  a 
generie  set  of  injeetor  element  types. 

2.0  Objectives 

The  overarching  objective  of  the  proposed  research  program  was  to  advance  our 
understanding  of  the  initiation  and  sustenance  of  eombustion  instability  in  rocket  engines  that 
use  a  liquid  oxygen/methane  propellant  eombination.  Utilizing  an  existing  multi-element  test 
chamber  that  is  described  in  the  next  sub-seetion,  we  propose  to  focus  on  the  following  set  of 
specific  objectives: 
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[1]  To  determine  the  ehamber  pressure  and  fuel  temperature  at  whieh  the  onset  of  transverse 
modes  due  to  eombustion  instability  are  observed  to  oeeur  for  suberitieal,  transeritieal 
and  supereritieal  propellant  eonditions  for  three  representative  injeetor  element  types. 

[2]  To  determine  a  pressure  sealing  relationship  for  eaeh  injeetor  element  for  the  onset  of 
eombustion  instability  based  on  measurements  done  over  a  range  of  diserete  ehamber 
pressures  and  propellant  temperature  eonditions. 

[3]  To  establish  a  baseline  set  of  data  to  eharaeterize  eombustion  instability  for  injeetors 
whieh  operate  over  a  10:1  throttling  range. 

[4]  To  transfer  the  understanding  gained  in  this  researeh  program  rapidly  to  industry  through 
elose  eooperation  with  the  leading  U.S.  liquid  propellant  roeket  manufaeturer,  Pratt  and 
Whitney. 

3.0  Experimental  Setup 

3.1  The  Cryogenic  Combustion  Laboratory 

A  series  of  eombustion  instability  studies  were  eondueted  at  the  Cryogenie  Combustion 
Laboratory  (CCL).  Several  publieations  eontain  information  regarding  the  layout  and 
eapabilities  of  the  CCL  [4-8].  Thus,  only  an  overview  of  the  faeilities  will  be  diseussed  here. 
The  CCL  faeility  eonsists  of  three  rooms  eonstrueted  with  reinforeed  eonerete  to  proteet 
operators  and  adjaeent  areas  from  any  hazardous  explosions  or  eonditions  whieh  may  oeeur. 
The  main  room  is  the  eontrol  room,  whieh  holds  the  data  aequisition  eomputers  and  flow  eontrol 
equipment.  A  programmable  logic  controller  (PLC)  remotely  aetuates  valves  and  monitors 
critieal  pressure  transdueers  in  the  flow  system.  If  a  pressure  transducer  or  other  sensor  reads 
outside  of  programmed  ranges,  an  abort  monitor  automatically  aborts  the  test  and  brings  on 
purge  gases  to  purge  the  system  of  combustible  propellants.  The  control  room  also  houses 
various  PC  computers  and  timing  systems  used  for  data  acquisition.  One  dedicated  eomputer 
located  at  the  control  panel  provides  the  operators  with  real  time  data  overviews  and  data 
acquisition  status.  Other  computers  are  used  for  other  acquisition  tasks,  so  that  critical 
performance  data  (pressure  and  temperature)  may  be  taken  separately  from  other  advanced 
diagnostics  (laser  based  techniques,  high  frequency  pressure  measurements,  etc.).  A  separate 
instrumentation  room  houses  lasers,  optics,  and  any  other  equipment  needed  for  data  capture 
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during  a  test  which  is  deemed  too  fragile  to  leave  next  to  the  rocket  should  a  hardware  failure 
occur.  Optical  and  cable  access  ports  are  available  to  route  the  necessary  equipment  or  sensors 
into  the  test  cell  for  data  capture. 

Flowrate  capabilities  at  the  CCL  provide  0.45  kg/s  (l.Olbm/s)  of  liquid  and  gaseous 
oxygen  (LOX  and  GOX),  0.11  kg/s  (0.25  Ibm/s)  of  gaseous  hydrogen  (GH2),  0.23  kg/s 
(0.50  Ibm/s)  of  liquid  hydrocarbon  (methanol,  RP-1,  etc.),  and  2.27  kg/s  (5.0  Ibm/s)  of  air. 
Propellant  flowrate  control  is  accomplished  using  critical  orifices  for  gaseous  flow  and  cavitating 
venturis  for  liquid  flows. 

3.2  Experimental  Hardware 

3.2. 1  Main  Chamber  Design 

Experiments  for  studies  discussed  in  this  report  were  conducted  in  a  multi-clement, 
rectangular  rocket  chamber.  Propellants  chosen  for  these  studies  were  liquid  oxygen  (LOX)  and 
gaseous  methane  (CH4).  Figure  1  shows  a  schematic  of  the  rocket  chamber  on  the  test  stand 
whereas  Fig.  2  shows  an  overhead  photograph  of  the  rocket  on  the  test  stand.  Flow  is  from  top 
to  bottom  in  the  picture.  When  assembled,  the  combustion  chamber  measures  203  mm  (8.00  in.) 
wide  by  81.3  mm  (3.20  in.)  long  (in  the  flow  direction)  by  19  mm  (0.75  in.)  high.  These 


Fig.  1.  Schematic  of  rectangular  chamber  on  test  stand. 
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Fig.  2.  Photograph  of  rocket  firing  (2  injectors,  8  nozzles,  Pc=222  psia,  0/F=3).  Flow  is  from 
top  to  bottom  in  the  photograph. 

dimensions  were  chosen  to  isolate  the  first  and  second  transverse  modes  (IW  and  2W  modes)  of 
the  chamber  at  frequencies  low  enough  (-2000-8000  Hz)  to  be  readily  studied.  At  the  head  end 
of  the  chamber  is  located  an  injector  body  with  ports  for  up  to  16  injector  elements.  At  the 
exhaust  end  of  the  chamber,  a  nozzle  body  is  located  with  ports  for  up  to  8  nozzle  elements. 
Each  injector  element  and  nozzle  element  could  be  independently  moved  within  the  injector 
body  and  nozzle  body,  respectively.  Injeetor/nozzle  element  refers  to  a  single  injeetor/nozzle, 
whereas  injeetor/nozzle  body  (or  block)  refers  to  the  component  that  holds  several 
injeetor/nozzle  elements  or  their  blanks,  and  is  multi-element  in  nature.  For  discussion 
throughout  the  report,  the  injectors  and  nozzles  have  been  numbered  for  clarity.  Figure  3  shows 
a  section  view  of  the  chamber  with  injector  and  nozzle  elements  numbered  and  the  igniter  ports 
and  water  passage  ports  identified 

3.2.2  Injector  and  Injector  Block 

The  primary  injector  clement  design  utilized  in  the  majority  of  studies  was  a  swirl  coaxial 
injector.  The  swirl  coaxial  injector  element  consisted  of  a  central  post  through  which  LOX  was 
introduced.  At  one  end  of  the  post  was  a  swirl  nut.  The  ID  of  the  post  was  2.36  mm  (0.093  in.) 
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Igniter  1 


Water 
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FLOW 
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Igniter  2 


Water 

Inlet 


Water-cooled  Nozzles 

Fig.  3.  Section  view  of  chamber  with  injector  elements  and  nozzle  elements  numbered. 


while  its  OD  was  3.30  mm  (0.130  in.).  The  post  length  was  50.8  mm  (2.00  in.)  from  the  end  of 
the  swirl  nut  to  the  injection  plane.  Six  tangential  holes  0.79  mm  (1/32  in.)  in  diameter  were 
located  in  the  swirl  nut.  The  holes  were  biased  three  to  a  side,  mirrored  about  the  central  axis  of 
the  post,  in  order  to  introduce  a  radial  component  to  the  flow.  The  injector  post  was  silver 
soldered  into  a  custom  made  fitting  which  was  used  to  insert  the  injector  clement  into  the  injector 
body.  The  fuel  annulus  was  formed  by  an  insert  on  the  injection  side  of  the  injector  body. 
The  primary  injector  clement  annulus  insert  was  3.91  mm  (0.154  in.)  outer  diameter.  Centering 
rings  located  between  the  annulus  insert  and  injector  body  kept  the  LOX  post  aligned. 
Stand-offs  of  various  lengths  could  be  inserted  into  the  upstream  side  of  the  injector  body  to 
produce  various  LOX  post  recesses.  For  the  experiments  discussed  here,  all  stand-offs  were 
5.08  mm  (0.200  in.)  in  length  to  produce  a  flush  injector  face.  A  schematic  of  the  injector 
clement  is  shown  in  Fig.  4  and  a  cut-away  schematic  showing  the  assembly  of  four  injector 
elements  in  the  injector  body  is  shown  in  Fig.  5. 

3.2.3  Propellant  Manifold 

To  ensure  single  phase  liquid  flow,  the  LOX  feed  lines  were  jacketed  and  cooled  with 
LN2.  Prior  to  connecting  with  the  injector  block,  the  LOX  manifold  is  submerged  in  a  LN2  bath 
with  independent  feed  lines  from  the  manifold  connecting  to  injector  elements. 
Each  independent  feed  line  also  has  a  LN2  jacket  to  keep  propellants  at  cryogenic  temperatures 
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Fig.  4.  Cut-away  schematic  of  injector  element. 


until  as  close  to  the  injection  block  as  possible.  After  the  metering  orifice,  a  large  manifold 
block  is  used  for  the  methane  propellant  with  independent  feed  lines  to  the  injector  block.  Ports 
for  the  methane  alternate  top  and  bottom  in  the  injector  block. 

Jensen,  Dodson  and  Claflin  [9,  10]  noted  that  temperature  ramping  was  an  effective 
means  of  rating  the  stability  characteristics  of  LOX/mcthanc  combustion.  Thus,  as  a  means  of 
causing  instabilities  in  the  chamber,  the  methane  was  reduced  temperature.  Special  hardware 
was  required  to  accomplish  this  scries  of  testing. 

In  order  to  cool  the  methane  gas  down  to  the  low  temperatures  required  for  testing  (-70  to 
-90°C),  a  special  bath  was  designed  for  the  methane  circuit.  The  bath  was  a  cryogenic  tank  filled 
with  a  specially  formulated  heat  transfer  fluid  from  Duratherm  Corp.  The  fluid  was  fomiulatcd 
to  remain  in  a  low  viscosity  liquid  state  with  high  heat  transfer  characteristics  even  at  reduced 
temperature.  An  emersion  cooler  (Haakc  EK-90)  was  then  used  to  bring  the  heat  transfer  fluid 
down  to  -90°C.  Propellant  lines  from  the  methane  manifold  were  then  coiled  through  the  bath 
to  create  a  simplified  heat  exchanger.  Additional  jackets  of  LN2  maintained  cooled  methane 
temperatures  to  the  injector  block. 

For  cooled  methane  testing,  the  run  valves  for  the  methane  circuit  were  located  upstream 
of  the  methane  manifold  block  and  feed  lines  from  the  manifold  were  independently  routed 


Fig.  5.  Cut-view  schematic  of  injector  block  showing  4  injector  elements. 
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through  the  methane  temperature  bath  to  the  injeetor  bloek.  Due  to  the  inereased  length  of  line 
for  the  methane  eireuit,  timing  for  the  valves  was  adjusted  to  aeeount  for  neeessary  line  fill  and 
pressurization.  Temperature  probes  in  the  methane  temperature  bath  and  on  two  feed  lines  just  at 
the  feed  line  exit  from  the  bath  provided  monitoring  of  methane  temperature.  A  LN2  eoil  from 
the  methane  temperature  bath  to  the  injeetor  bloek  provided  a  triekle  flow  of  liquid  nitrogen  to 
overeome  temperature  loss  from  the  bath  to  the  injeetor  bloek. 

3.2.4  Nozzle  and  Nozzle  Bloek 

Nozzle  elements  were  all  of  a  singular  design,  fabrieated  from  molybdenum,  whieh  eould 
be  mounted  in  a  nozzle  bloek  permitting  up  to  maximum  of  8  nozzles  elements.  Eaeh  nozzle 
element  was  eonieal  in  design  with  a  throat  diameter  of  5.44  mm  (0.214  in.).  When  the  nozzle 
element  was  inserted  into  the  nozzle  bloek,  a  flow  passage  was  formed  around  the  throat  region 
nozzle  element  that  was  eooled  using  high  pressure  water  as  the  eoolant.  When  fewer  than  eight 
nozzle  elements  were  being  used,  blanks  with  the  same  shape,  but  no  throat  eut-out,  were 
inserted  in  empty  spaees. 

3.2.5  Igniters 

Two  GO2/GH2  igniters  loeated  on  the  side  walls  of  the  ehamber  provided  uniform 
ignition  aeross  the  ehamber,  and  are  the  same  as  those  diseussed  in  other  experiments  [4,  5]. 
The  igniters  were  fed  from  gas  bottles  independent  of  the  main  propellant  feed  system. 

3.2.6  High  Frequenev  Pressure  Transdueers 

PCB  Eleetronies  model  1 13A24  high  frequeney  pressure  transdueers  were  used  for  high 
frequeney  pressure  measurements  in  the  ehamber.  Two  sensors  were  used  in  the  ehamber,  one 
mounted  along  the  eenterline  of  the  ehamber,  the  other  near  the  side  wall.  Positions  were  ehosen 
to  measure  the  pressure  oseillations  the  anti-nodes  of  the  IW  and  2W  modes  (for  the  ease  of  the 
wall  pressure  transdueer),  as  well  as  the  node  of  the  IW  mode  and  anti-node  of  the  2W  mode 
(for  the  ease  of  the  eenterline  pressure  transdueer).  Figure  6  shows  the  positions  of  the  sensors 
with  respeet  to  the  mode  shapes  of  the  IW  and  2W  modes.  These  sensors  have  a  resonant 
frequeney  of  500  kHz  and  a  pressure  range  of  0-1000  psi.  Samples  were  taken  at  100  kHz,  so  a 
maximum  useful  frequeney  range  of  0-50  kHz  was  available.  The  transdueers  were  mounted  in 
water  eooled  jaekets  and  reeessed  into  the  ehamber  wall  to  proteet  them  from  the  hot  eombustion 
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Fig.  6.  Schematic  of  sensor  positions  with  respect  to  mode  shapes. 

gases.  Openings  to  the  chamber  were  1.59  mm  (0.063  in.)  diameter,  12.2  mm  (0.480  in.)  long 
for  the  centerline  transducer,  and  18.5  mm  (0.730  in.)  long  for  the  wall  transducer.  Difference  in 
measurement  port  lengths  was  due  to  the  necessity  of  keeping  the  wall  transducer  recessed  back 
further  so  the  port  opening  could  be  as  close  to  the  side  wall  of  the  chamber  as  possible. 
Recesses  were  kept  as  much  to  a  minimum  as  hardware  allowed,  in  order  to  prevent  spurious 
modes  from  being  detected. 

3.3  Test  Conditions 

3.3. 1  Basic  Flow  Conditions 

Each  injector  clement  was  designed  based  on  a  LOX  flowrate  of  0.057  kg/s  (0. 125  Ibm/s) 
per  clement.  The  overall  mixture  ratio  had  an  0/F  =  3  ((|)  =  1 .33),  with  methane  (CH4)  as  the 
fuel.  A  total  flowrate  capability  of  0.454  kg/s  (1 .00  Ibm/s)  of  LOX  flow  was  used. 
For  reference,  LOX/methane  stoichiometric  mixture  ratio  is  0/F  =  4.  Table  1  lists  the  critical 
pressure  and  temperatures  for  LOX  and  methane  for  reference  in  evaluating  the  thermodynamic 
state  of  the  propellants  for  the  current  studies.  Nozzle  elements  could  be  altered  to  change  the 
chamber  pressure  to  the  desired  quantity. 


Table  1.  Table  of  critical  temperature  and  pressure  for  LOX  and  methane. 


Propellant 

Critical  Temp.  (K  [°F1) 

Critical  Pressure  (MPa  [psia]) 

Liquid  Oxygen  (LOX) 

154.6  [-181.4] 

5.043  [731.4] 

Methane 

190.6  [-116.6] 

4.599  [667.1] 
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Fig.  7.  Sample  pressure  trace  for  mean  chamber  pressure,  indicating  the  steady  state  period  used 
during  test  runs. 

3.3.2  Basic  Timing  Sequence 

The  testing  sequence  for  the  rocket  chamber  consists  of  two  stages,  a  start-stage,  and  a 
run-stage.  The  start  stage  was  included  to  prove  the  ignition  of  the  propellant  flows  as  well  as  to 
case  the  pressure  shock  to  the  system  as  it  comes  up  to  full  pressure.  Since  methane  is  more 
difficult  to  successfully  ignite,  the  start  stage  utilized  hydrogen  to  ensure  ignition,  and  then  the 
hydrogen  was  shut  off  as  methane  was  brought  online  during  the  transition  to  the  run  stage. 
Figure  7  shows  a  sample  pressure  trace  for  a  typical  test.  Pressures  are  shown  in  psia,  and  are 
referenced  to  test  run  time. 

4.0  Results  and  Discussion 

4. 1  High  Frequency  Data  Analysis  Overview 

4. 1 . 1  Chamber  Aeoustie  Characteristies 

Before  presenting  the  results,  an  overview  of  the  analysis  method  for  the  high  frequency 
pressure  transducer  measurements  is  presented.  The  modal  calculations  for  the  chamber  are 
based  upon  solution  of  the  wave  equation  for  a  rectangular  cavity.  The  equation  describing  the 
aeoustie  modal  frequencies  of  the  cavity  is  given  by: 
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Fig.  8.  Sehematies  illustrating  coordinate  system  in  describing  acoustic  modes  of  chamber. 
Left  image:  Section  view  of  chamber  with  included  dimensions.  Right  image;  Schematic  of 
chamber  on  test  stand  with  included  coordinate  system. 
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where  a  is  the  acoustic  velocity  in  the  chamber,  L^,  Ly,  and  L~  are  the  dimensions  of  the  chamber 
(see  Fig.  8)  and  /,  m  and  n  are  the  mode  numbers  corresponding  to  the  jc,  y,  and  z  coordinates, 
respectively.  In  this  analysis,  /  corresponds  to  the  transverse  mode  related  to  Ly,  m  to  the 
longitudinal  mode  related  to  Ly,  and  n  corresponds  to  the  transverse  mode  related  to  L:. 
The  sehematies  in  Fig.  8  illustrate  the  dimensions  of  the  chamber  {Ly,  Ly,  and  L^)  with  respect  to 
the  orientation  of  the  chamber  on  the  test  stand.  The  left  side  image  is  a  section  view  of  the 
combustion  chamber  with  dimensions  included,  whereas  the  right  side  image  shows  the  rocket 
chamber  assembly  as  placed  on  the  test  stand  with  the  coordinate  system  shown.  In  these 
images,  Ly  is  the  width  of  the  chamber  (203  mm  [8  in.]),  Ly  is  the  length  (81.3  mm  [3.2  in.]),  and 
L;  is  the  hcight(19.1  mm  [0.75  in.]).  In  referring  to  the  modes  of  the  chamber,  the  nomenclature 
for  the  modes  in  the  transverse,  longitudinal  (or  axial),  and  height  {2nd  transverse)  dimensions 
are  IW,  mL,  and  nH,  respectively. 


The  dimensions  of  the  rectangular  chamber  were  selected  such  that  the  frequencies  for  the 
longitudinal  and  2nd  transverse  modes  were  much  greater  than  the  IW  (transverse)  mode. 
Table  2  lists  the  expected  modal  frequencies  of  the  chamber  assuming  an  equilibrium  acoustic 
velocity  (calculated  from  CEA  [11,  12])  of  1277  m/s  (4188  ft/s).  CEA  calculations  for  the 
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Table  2.  Tabic  of  expected  mode  frequencies  (Hz). 


Mode  No. 

x=\ 

It 

jc  =  3 

xW 

3141 

6282 

9423 

xL 

7853 

15,706 

23,559 

xH 

33,507 

67,013 

100,521 

chamber  with  inputs  of  888  psi  chamber  pressure  and  0/F  =  3  were  used  to  determine  the 
chamber  gas  acoustic  velocity.  The  difference  in  acoustic  speed  for  222  psi  tests  is  only  2% 
lower  (1250  m/s),  and  thus  does  not  significantly  change  expected  mode  frequencies  (3078  Hz 
vs.  3141  Hz  for  IW  mode).  Note  that  the  frequencies  for  the  IL  and  IH  shown  in  Table  2  are 
well  above  the  IW  and  2W  modes  of  the  chamber.  Thus,  the  chamber  is  well  suited  to  isolating 
the  IW  and  2W  modes  for  the  present  study.  Given  that  fact  and  considering  that  the  transverse 
modes  are  of  primary  interest.  Equation  1  may  be  simplified  to  a  1-D  equation,  given  by: 


J 


[2] 


It  was  observed  during  testing  that  the  predicted  modal  frequency  response  did  not  match 
exactly  with  the  frequencies  given  in  Table  2.  However,  the  frequencies  were  consistently 
shifted  from  equilibrium  calculations  to  slightly  lower  values  (c.g.  2500-2700  Hz  as  compared  to 
the  CEA  predicted  3141  Hz).  While  the  exact  cause  of  this  difference  is  unknown,  departures 
from  a  uniform  gas  temperature  field  in  the  chamber  and  reduced  c*  efficiencies  (?/c*  <  1)  are 
suspected  as  reasons  for  reduced  gas  temperature  in  the  combustion  chamber  and  thus  the 
observed  lower  frequency  response. 

4.1.2  Power  Spectrum  Density  Overview 

Since  the  modal  response  of  the  chamber  is  frequency  dependent,  a  method  is  needed  to 
analyze  from  the  high  frequency  signal  measured  in  the  time  domain  and  transform  it  to  the 
frequency  domain.  In  the  current  study,  the  high  frequency  pressure  signal  is  transformed  into  a 
Power  Spectrum  Density  (PSD)  plot.  The  PSD  is  related  to  the  Fourier  transform  as  shall  be 
discussed  next. 


The  most  common  and  readily  available  method  for  this  frequency  analysis  is  the  fast 
Fourier  transform  (FFT).  The  Fourier  transform  (FT)  is  given  by: 
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where  J{t)  is  the  value  of  the  fiinetion  (or  signal)  in  the  time  domain  where  t  is  time,  /  is  the 
imaginary  number  ( ),  co  is  the  angular  frequeney  ( o)  =  l7f ),  and  /  is  the  value  of  the 
funetion  (or  signal)  transformed  into  the  frequeney  domain. 

The  fast  Fourier  transform  is  the  diseretized  form  of  the  FT,  written  to  take  advantage  of  a 
eomputer's  eomputational  abilities  and  the  nature  of  a  diseretized  sample.  Its  form  is  given  by: 

)t  =  0,l,2,...A^-l  [4] 

n*0 

Here,  X  and  x  are  the  values  of  the  funetion  in  the  frequeney  and  time  domains,  respeetively.  N  is 
the  total  number  of  data  points  to  be  sampled  over  some  time  interval,  where  N  is  restrieted  to 
powers  of  2.  A/  is  the  inerement  of  time,  or  sampling  period  ( 1 /sampling  rate).  Both  k  and  n  are 
indiees  used  in  the  summation.  A/is  the  frequeney  resolution,  defined  as  the  sampling  frequeney 
divided  by  the  total  number  of  samples,  N. 

Related  to  the  Fourier  transform,  the  power  speetrum  density  (PSD)  is  a  measure  of  how 
the  power  of  a  signal  is  distributed  in  the  frequeney  domain.  The  power  speetrum  of  a  time 
dependent  signal, y(x),  is  given  by: 

5x.(/)=^*(/k{/)  =  |4/f  [5] 

where  X(J)  is  the  Fourier  transform  ofy(jc)  and  X*if)  is  the  eomplex  eonjugate  of  X{J).  Thus  a 
plot  of  the  power  speetrum  with  respeet  to  frequeney  is  termed  the  power  speetrum  density. 
Integrating  the  power  speetrum  over  the  entire  frequeney  domain  gives  the  total  intensity 
(expressed  as  psi^Hz),  or  total  varianee  (a^,or),  of  the  entire  signal,  as  expressed  in  Eq.  6. 

^l.  =  ^f{copco  [6] 

Integrating  Eq.  (4.6)  from  o  to  6  as  opposed  to  the  entire  domain  gives  the  intensity  of  the  signal 
between  the  frequeneies  of  a  and  h. 

4.1.3  Power  Speetrum  Density  Analysis 

Beeause  of  the  properties  of  Eq.  6,  the  PSD  provides  an  effeetive  means  of  analyzing  the 
high  frequeney  data.  Integrating  the  PSD  over  the  frequeney  range  between  f\  to  72  gives  the 
intensity,  also  ealled  the  varianee  (a  ),  of  the  signal  within  that  frequeney  band.  As  noted  in  the 
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previous  section  (4. 1 .2),  this  is  related  to  the  pressure  oscillation  amplitude  within  this  frequency 
band.  Taking  the  square  root  of  the  variance  then  gives  the  root  mean  squared  amplitude  of  the 
pressure  oscillations  (P'rms)  within  the  desired  frequency  band.  Multiplying  the  rms  pressure 
oscillation  amplitude  by  2^2  yields  the  pcak-to-peak  amplitude  of  the  pressure  oscillations 
(P'p.p).  It  is  then  possible  to  relate  the  amplitude  of  pressure  oscillations  within  a  given  frequency 
band  relative  to  mean  chamber  pressure  (P'^yPe). 

A  set  of  virtual  instrument  (VI)  written  for  the  LabVlEW  software  provides  a  method  for 
quickly  obtaining  the  PSD  of  the  high  frequency  pressure  signal.  By  measuring  the  rms  and 
pcak-to-peak  amplitudes  of  modal  frequencies  in  the  signal,  and  normalizing  by  mean  chamber 
pressure,  the  mean  instability  level  (defined  here  as  Pp.^/Pc  and  expressed  as  a  percentage)  can  be 
determined.  This  level  is  expressed  as  a  percentage  of  the  mean  chamber  pressure  for  a  given 
mode  over  any  time  interval  of  interest.  The  program  can  also  be  used  to  obtain  an  array  of  PSD 
variance  and  peak-to-pcak  pressure  amplitude  for  a  given  frequency  band  as  a  function  of  time. 
Thus  a  plot  of  the  peak  to  peak  pressure  amplitude  versus  time  can  be  obtained.  By  comparing 
plots  from  different  tests  where  conditions  such  as  injector  clement  placement,  flow  conditions 
and  nozzle  clement  placement  were  changed,  stability  characteristics  may  be  obtained  based  on 
these  changes  in  test  conditions. 

4.2  Chamber  Response 

4.2.1  Acoustic  Response  Studies 

The  initial  set  of  experiments  conducted  in  the  rectangular  rocket  chamber  focused  on 
measuring  the  chamber  frequency  response  as  a  function  of  injector  placement  and  chamber 
pressure.  This  scries  of  experiments  measured  the  acoustic  response  to  determine  whether 
spontaneous  combustion  instability  could  be  observed,  as  no  external  methods  were  used  to 
excite  an  instability  level  in  the  chamber.  By  placing  injectors  at  pressure  anti-nodes,  it  was 
expected  that  instabilities  could  be  readily  generated  as  compared  to  those  tests  where  injectors 
were  placed  at  pressure  nodes.  It  could  then  be  observed  whether  spontaneous  instabilities 
occurred  within  the  chamber  or  not. 

As  discussed  in  See.  4.1.2,  the  primary  method  of  analysis  of  the  high  frequency  pressure 
signal  was  the  use  of  the  PSD.  Through  the  PSD,  it  was  possible  to  identify  which  modes  were 
active  in  the  chamber,  and  the  relative  amplitude  of  those  signals.  Figure  9  shows  a  sample  PSD 
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Fig.  9.  PSD  plot  of  wall  transducer  for  2  injector  elements  at  the  walls  (Pos.  1  and  16), 
Pc  =  222  psia. 

plot  of  the  pressure  transducer  located  near  the  wall  (ref.  Fig.  6  for  position  of  the  pressure 
transducer)  for  the  ease  of  two  injectors  located  at  positions  I  and  16  (ref  Fig.  3)  and  a  target 
chamber  pressure  of  1.53  MPa  (222  psi).  The  PSD  shown  was  taken  for  2'^  data  points  (16384 
data  points  or  0.16  s  of  steady  state  portion  of  the  pressure  signal).  In  this  figure  the  x-axis  is 
frequency,  whereas  the  amplitude  (y-axis)  indicates  the  PSD  amplitude  (psi^Hz).  Figures  9  (for 
the  wall  pressure  transducer)  and  1 0  (for  the  centerline  pressure  transducer)  both  have  the  1 W 
and  2W  modes  at  approximately  2500  Hz  and  5000  Hz,  respectively.  For  comparison,  the  ideal 
modal  frequencies  based  on  equilibrium  conditions  are  3 101  Hz  and  6022  Hz  for  the  1 W  and  2W 
modes,  respectively.  Departures  from  these  ideal  frequencies  are  likely  due  to  differences  in  gas 
temperature  within  the  chamber,  as  evidenced  by  the  approximate  c*  efficiencies  which  ranged 
from  70-102%  for  all  tests  conducted.  The  mean  and  median  values  of  c*  for  all  tests  conducted 
were  both  91%.  It  is  also  clear  that  the  1 W  mode  PSD  amplitude  is  quite  strong  compared  to  the 
2W  mode  PSD  amplitude,  indicating  that  most  of  the  acoustic  energy  is  present  in  the  1 W  mode. 

By  comparison.  Fig.  10  shows  an  example  of  the  PSD  for  the  same  test  runs,  for  the 
centerline  pressure  transducer.  The  axes  of  the  plot  arc  the  same  as  previously  described,  and  the 
data  sample  size  is  the  same.  For  this  ease,  it  is  noted  that  the  1 W  mode  is  significantly  reduced 
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Fig.  10.  PSD  plot  of  centerline  transducer  for  2  injector  elements  at  the  walls  (Pos.  1  and  16), 
Pc  =  222  psia. 

compared  to  the  wall  pressure  transducer  plot  (Fig.  9).  This  is  expected,  as  the  IW  mode  exhibits 
a  pressure  node  at  the  centerline  position. 

From  observation,  it  is  noted  that  the  amplitude  of  the  2W  mode  is  not  the  same  for  the 
wall  and  centerline  transducers,  despite  the  symmetric  nature  of  the  2W  mode.  Ideally,  these 
amplitudes  should  be  the  same,  and  no  source  for  this  disparity  has  been  identified.  One  possible 
explanation,  however,  is  due  to  local  differences  in  the  gas  temperature  and  composition  causing 
local  differences  in  acoustic  properties  and  modal  eharaeteristies,  although  it  is  not  possible  from 
the  data  to  confirm  this.  Based  on  static  pressure  transducer  measurements,  the  actual  chamber 
pressure  was  lower  than  the  target  condition,  with  an  approximate  c*  efficiency  for  this  test 
condition  of  98%.  It  is  worth  noting  here  that  acoustic  characteristics  of  the  chamber  are 
assumed  on  linear  acoustic  theory  in  a  confined  duet  with  no  additional  openings. 

4.2.2  Burst  Phenomenon 

Upon  further  examination  of  the  raw  high  frequency  pressure  signals,  it  was  observed  that 
pressure  oscillations  in  the  chamber  grew  and  decayed  in  amplitude  quickly.  These  oscillations 
occur  in  random  patterns  throughout  the  test  sequence  and  are  hereafter  referred  to  as  “burst" 
phenomenon.  Figure  1 1  illustrates  the  burst  phenomenon  within  the  chamber.  The  main  plot  of 
the  figure  shows  a  sample  time  trace,  in  this  ease  for  a  two  injector  ease  (Pos.  1,16),  four  nozzles. 
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Fig.  11.  Illustration  of  burst  phenomenon  during  test  series. 


target  P^.  =  444  psia,  with  0. 143  in.  l.D.  annuluses  (higher  momentum  ratio).  The  y-axis  shows 
the  amplitude  of  the  pressure  signal  above  mean  ehamber  pressure,  whereas  the  x-axis  shows  the 
data  point  number  of  the  cycle.  The  x-axis  is  also  representative  of  the  time  scale,  since  the 
sample  rate  was  1x10^  pts./s.  To  obtain  the  time  scale,  simply  divide  the  point  number  by  the 
sample  rate.  Thus,  while  the  x-axis  begins  at  340000,  this  is  analogous  to  3.4  s.  The  main  plot 
of  the  figure  represents  0.1  s  worth  of  data  from  the  steady  state  portion  of  the  rocket  firing. 
At  this  scale,  the  burst  phenomenon  of  the  instability  is  evident.  The  inset  of  the  figure  is  a 
zoomed  in  portion  of  the  main  graph,  showing  a  segment  approximately  0.02  s  in  duration. 
The  IW  mode  frequency  is  clearly  evident  in  this  inset.  Examination  of  the  figure  shows  that 
IW  mode  oscillation  amplitudes  (P'p.;,)  during  this  segment  grow  from  approximately  10  psi 
(~  2.5%  P'p.p/Pc),  to  a  level  of  40  psi  (~10%  P'p.p/Pc).  The  pattern  of  the  bursts  is  random  in 
nature.  In  the  inset,  the  main  burst  shown  is  approximately  12  cycles  in  duration,  with  a 
predominant  1 W  mode  frequency  of -2500  Hz. 

4.2.3  Short-Time  Fourier  Analysis 

Because  of  the  burst  nature  of  the  instabilities  in  the  ehamber,  the  previous  method  of 
examining  instabilities  in  the  ehamber  as  an  average  segment  (ref  See.  4.2.1)  was  not  sufficient 
to  fully  describe  the  instability  level  in  the  ehamber.  Thus,  a  short-time  Fourier  analysis  (STFA) 
was  performed  to  provide  a  time  history  of  the  instability  level  in  the  chamber.  In  this  method,  a 
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Fig.  12.  Example  plot  of  1 W  mode  vs.  time. 


FFT  was  run  on  the  signal  for  a  short  period  of  time  (in  this  ease  512  data  points),  then  the  FFT 
analysis  advaneed  forward  by  500  data  points,  providing  a  small  overlap  between  subsequent 
FFTs  to  a  new  segment  of  signal.  Throughout  the  analysis,  the  varianee  and  peak-to-peak 
oseillation  amplitude  is  measured  for  a  partieular  mode  of  interest  (either  the  1 W  or  2W  mode). 
Thus,  a  running  value  of  the  varianee  and  instability  level  of  partieular  modes  as  a  funetion  of 
time  ean  be  obtained.  Figure  12  shows  an  example  of  a  IW  mode  plot  obtained  by  STFA.  In 
Fig.  12,  the  y-axis  is  the  instability  level  (P'p.p/?c)  and  the  x-axis  is  the  runtime  (note  steady  state 
ehamber  pressure  data  oeeurs  from  3.2  to  4.0  seeonds  during  the  test  sequenee).  Several  runs  at 
the  same  test  eonditions  are  shown.  The  blue  dashed  line  in  the  figure  represents  the  10% 
instability  level,  whereas  the  red  dashed  line  represents  the  5%  instability  level,  shown 
for  quiek  referenee.  Blaek  dashed  lines  braeket  the  approximate  maximum  and  minimum 
instability  levels  of  the  partieular  test  eondition.  Thus,  STFA  plots  for  various  modes  may  be 
quiekly  eompared  to  see  how  ehanges  in  test  eonditions  (sueh  as  injeetor  plaeement)  affeet 
instability  levels. 
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4.2.4  Temperature  Ramping  Studies 


Since  the  pulsation  device  did  not  induce  instability  levels  (P'pyPf)  of  the  magnitude 
expected,  other  methods  of  driving  instabilities  were  investigated.  In  partieular,  following  the 
work  of  Jensen,  et  al.  [9,  10],  reduetion  of  the  methane  gas  temperature  appeared  a  possible 
mechanism  for  inducing  instabilities  within  the  combustion  chamber  and  was  of  particular 
interest  for  the  this  study.  The  method  utilized  in  the  present  work  for  reducing  methane  gas 
temperature  used  a  simple  thermal  heat  exchanger.  Based  on  the  studies  by  Jensen,  et  al.,  [9],  it 
was  determined  that  a  target  methane  gas  temperature  of  -90'’C  would  be  sufficient  for  inducing 
instabilities.  Tests  were  conducted  for  both  ambient  (room  temperature)  and  cooled  methane 
runs  for  the  same  chamber  conditions.  Figure  13  shows  a  sample  IW  mode  plot  for  a  cooled 
methane  test,  whereas  Fig.  14  shows  a  IW  mode  plot  for  an  ambient  temperature  test. 
The  y-axis  and  x-axis  are  the  same  as  those  discussed  for  Fig.  12.  The  tests  shown  were  for  four 
injector  elements  (Pos.  13,  14,  15,  16)  with  eight  nozzles  (P^.  =  444  psia),  as  illustrated  by  the 
graphic  in  the  upper  left  comer.  It  should  be  noted  that  the  mean  methane  gas  temperature, 
measured  just  upstream  of  injection  was  -70  to  -100°F  (-56  to  -73°C)  for  the  cooled  studies  and 
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Fig.  13.  Example  plot  of  IW  mode  vs.  time  for  cooled  methane  tests. 
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Fig.  14.  Example  plot  of  IW  mode  vs.  time  for  ambient  temperature  methane  tests. 


was  -6  to  +40°F  (-21  to  +4®C)  for  the  ambient  temperature  studies.  From  these  results,  one 
major  eonelusion  ean  be  drawn.  Despite  evidenee  from  Jensen,  et  ai  [9]  that  eooled  methane  is 
a  viable  option  for  performing  stability  rating  of  a  ehamber,  it  was  insuffieient  in  the  eurrent 
studies  to  show  any  significant  effect  on  the  chamber.  In  fact,  a  few  tests  conducted  with  cooled 
methane  exhibited  even  more  stable  operation  than  those  with  ambient  temperature  propellants 
(as  shown  in  Figs.  13  and  14).  However,  it  is  unclear  as  to  why  this  current  study  shows  results 
in  contradiction  with  those  of  Jensen,  et  al  [9].  However,  there  arc  some  key  differences 
between  the  tests.  One,  Jensen,  et  al  [9]  note  that  the  instabilities  observed  during  the  cooled 
tests  were  not  necessarily  IW  modes,  but  rather  an  odd  mode  that  did  not  correspond  to  chamber 
acoustics.  Finally,  the  work  by  Jensen,  et  al.  [9]  utilized  shear  coaxial  injectors  while  the  current 
work  used  swirl  coaxial  injectors.  While  similar  in  concept,  there  are  inherent  differences 
between  the  two  injector  types  including  mixing  and  atomization  properties.  Jensen,  et  al.  [9] 
note  that  the  odd  mode  may  be  related  to  organ  pipe  modes  of  the  fuel  annulus  with  the 
combustion  chamber.  Thus,  a  meaningful  one-to-one  comparison  between  the  current  and 
previous  bodies  of  work  may  not  be  appropriate. 
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Fig.  15.  Example  plot  of  IW  mode  vs.  time  for  cooled  methane  tests-showing  different 
instability  levels  for  similar  test  conditions. 


One  of  the  most  common  problems  with  combustion  instability  is  the  elusive  and 
unpredictable  nature  of  occurrences.  Simply  having  conditions  which  favor  occurrences  of 
instabilities  does  not  guarantee  that  an  event  will  occur.  Thus,  as  noted  by  Sutton  [13],  it  may 
take  numerous  tests  to  detect  a  potentially  catastrophic  condition.  The  current  scries  is  no 
exception  to  this  trend.  While  most  test  conditions  conducted  in  the  current  studied  showed 
similar  behavior  run  to  run,  a  few  test  conditions  showed  dramatically  different  behavior  from 
run  to  run.  Figure  15  shows  an  example  of  a  test  condition  whose  individual  runs  exhibited 
dramatically  different  results.  The  figure,  which  is  a  plot  of  the  IW  mode  similar  to  those 
discussed  previously,  shows  5  test  runs  for  the  chamber  operating  with  two  injector  elements 
located  at  the  wall  (Pos.  14  and  16),  eight  nozzles,  a  target  chamber  pressure  of  ?£.  =  222  psia, 
with  cooled  methane  propellant  (T,„y  = -105°F).  Of  the  five  runs  shown,  two  had  rather 
significant  instability  levels  greater  than  5%,  two  had  moderate  levels  between  3-5%,  and  one 
had  negligible  instability  levels. 
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Table  4.3.  Summary  of  test  conditions  showing  differing  instability  levels  in  chamber. 


Run 

# 

Chamber  Pressure 
MPa  Ipsia] 

Methane  T4 
°C  [°F1 

LOX  T4 
°C  |°F1 

O/F 

Instability  Level  (%) 
(P’pVPc) 

4 

1.43  [207.9] 

-73.5  [-100.3] 

-169.9  [-273.8] 

2.93 

4.22 

5 

1.45  [210.1] 

-75.6  [-104.0] 

-172.3  [-278.1] 

2.99 

6.76 

6 

1.45  [210.4] 

-76.2  [-105.1] 

-173.9  [-281.0] 

3.04 

6.80 

7 

1.46  [211.1] 

-76.7  [-106.0] 

-174.4  [-282.0] 

3.04 

3.97 

8 

1.45  [210.7] 

-76.9  [-106.4] 

-174.7  [-282.5] 

3.04 

2.05 

The  tests  were  all  conducted  on  the  same  day,  under  similar  flow  conditions,  with  no 
significant  departures  from  target  flow  conditions.  Actual  chamber  pressure  for  these  tests  were 
approximately  210psia  on  average,  with  an  approximate  c*  effieieney  of  90%.  The  tests,  in 
order,  had  instability  levels  of  4.22,  6.76,  6.80,  3.97,  and  2.05%  mean  instability  lcvcl(P'p.p/Pc). 
Table  3  lists  the  significant  flow  conditions  for  these  tests. 

5.0  Summary  and  Conclusions 

Studies  were  conducted  using  a  rectangular  rocket  chamber  burning  LOX/methane  to 
investigate  the  effect  of  methane  temperature  on  combustion  instability  initiation.  LOX  centered 
swirl  injector  elements  were  utilized  in  all  the  tests.  Chamber  pressures  between  222  and 
444  psia  were  studied  by  varying  the  number  of  injector  elements,  injector  placement  and 
methane  temperature.  Hot-fire  tests  with  injector  elements  positioned  near  the  pressure  anti- 
node  of  the  first  transverse  mode  exhibited  the  highest  pressure  oscillations. 

The  methane  propellant  temperature  was  varied  over  a  temperature  range  from  -70  to 
-100°F  (-56  to  -73°C)  for  the  cooled  studies  and  from  -6  to  +40°F  (-21  to  -l-4°C)  for  the  ambient 
temperature  studies.  The  results  obtained  from  monitoring  the  pressure  oscillations  using  high 
speed  pressure  transducers  indicated  no  difference  between  the  cooled  methane  test  eases  as 
compared  to  the  ambient  temperature  test  eases.  Pressure  oscillations  were  observed  for  some 
tests  to  be  as  high  as  10%  of  the  mean  chamber  pressure  for  both  the  cooled  and  ambient 
temperature  eases.  However,  repeated  tests  at  the  same  conditions  showed  the  magnitude  of  the 
observed  pressure  oscillations  varied  from  2-10%.  Thus,  the  high  variability  in  the  observed 
pressure  oscillations  appears  to  have  no  relationship  to  the  temperature  of  the  methane 
propellant. 
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